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f*"*. ' We discuss the calculation of the n/-contributions to the two-loop am- 

. plitude for e + e~ — * qgq. The calculation uses an efficient method based 

on nested sums. The result is presented in terms of multiple polyloga- 
rithms with simple arguments, which allow for analytic continuation in a 
■ straightforward manner. 
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<D 1. Introduction 

43 " 

Searches for new physics in particle physics rely to a large extend on our 
ability to constrain the parameters of the standard model. For instance, 
the strong coupling constant a s can be measured by using the data for 
e + e~ — ► 3-jets. At present, the error on the extraction of a s from this mea- 
surement is dominated by theoretical uncertainties [1], most prominently, 
by the truncation of the perturbative expansion at a fixed order. 

The perturbative QCD calculation of e + e~ — * 3-jets at next-to-next- 
to-leading order (NNLO) requires the tree-level amplitudes for e + e~ — ► 
5 partons [2], the one-loop amplitudes for e + e~ — ► 4 partons [3, 4] as well 
as the two-loop amplitude for e + e~ — ► qqg together with the one-loop ampli- 
tude e + e~ — > qqg to order e 2 in the parameter of dimensional regularization. 

The helicity averaged squared matrix elements at the two-loop level for 
e + e _ — ► qqg have recently been given [5]. In contrast, having the two-loop 
amplitude available, one keeps the full correlation between the incoming 
e + e~ and the outgoing parton's spins and momenta. Thus, one can study 
oriented event-shape observables. In addition, one has the option to inves- 
tigate event-shape observables in polarized e + e~-annihillation at a future 
linear e + e~-collider TESLA. 
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2. Calculation 

We are interested in the following reaction 

e + +e- ^q + g + q, (1) 

which we consider in the form, — > q(pi) + g{p2) + q~(P3) + e~{Pi) +e + {Pb), 
with all particles in the final state, to be consistent with earlier work [3]. 
The kinematical invariants for this reaction are denoted by 

Sij = {Pi + Pj f , S ijk = {pi + Pj + Pk) 2 , S = S123 , (2) 

and it is convenient to introduce the dimensionless quantities 

Xi = , X 2 = • (3) 

S123 $123 

Working in a helicity basis, it suffices to consider the pure photon ex- 
change amplitude A-y as it allows the reconstruction of the full amplitude 
with Z-boson exchange by adjusting the couplings. Furthermore, the com- 
plete information about A~, is given by just one independent helicity ampli- 
tude, which we take to be Ay(l + , 2 + , 3 _ , 4 + , 5~). All other helicity config- 
urations can be obtained from parity and charge conjugation. 

We can write yl 7 (l + , 2 + , 3 _ , 4 + , 5~) in terms of coefficients C2, C4, cq and 
c\2 for the various independent spinor structure as 

^ 7 (l + ,2+,3-,4 + ,5-) = -^M (4) 

(1 - Xi) [c 2 + ^C 6 - Ci 2 ^ + (1 - X 2 ) (c 4 - C12) + 2ci2 



x |s(35)[42] 
-<31>[12] 



[43] (35) (^2 + ^c 6 - cia) + [41] (15) (c 4 - c 12 ) J , 

where we have introduced the short-hand notation for spinors of definite 
helicity, \i±) = \pi±) = u±(pi) = v T (pi), (i ± | = (pi ± \ = u±(pi) = v T {p-i), 
and for the spinor products {pq) = (p — |g+) and [pq] = (p + \q— ). 

The coefficients q depend on the x\ and X2 of eq.(3) and can be calcu- 
lated in conventional dimensional regularization. To that end, we proceed 
as follows [6] -[8]. In a first step, with the help of Schwinger parameters [9], 
we map the tensor integrals to combinations of scalar integrals in various 
dimensions and with various powers vi of the propagators. For every basic 
topology, these scalar integrals can be written as nested sums involving T- 
functions. The evaluation of the nested sums proceeds systematically with 
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the help of the algorithms of [6] , which rely on the algebraic properties of 
the so called Z-sums, 
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Z(n;m 1 ,...,m k ;x 1 ,...,x k ) = ^ — ^ . . . . (5) 

n>ii>i 2 >...>jfc>0 1 k 

By means of recursion the algorithms allow to solve the nested sums in 
terms of a given basis in Z-sums to any order in e. .Z-sums can be viewed 
as generalizations of harmonic sums [10] and an important subset of .Z-sums 
are multiple polylogarithms [11], 

Lim fc ,...,miOEjfc, ...,xi) = Z(oo;mi,...,m k ;xi,...,x k ). (6) 

All algorithms for this procedure have been implemented in FORM [12] 
and in the GiNaC framework [13, 14]. In this way, we could calculate all 
loop integrals contributing to the one- and two-loop virtual amplitudes very 
efficiently in terms of multiple polylogarithms. 

The perturbative expansion in a s of the functions q is defined through 

c, = ^(<f + (£)<£> +(£)>+<>(«.•)) . (7) 

Then, after ultraviolet renormalization, the infrared pole structure of the 
renormalized coefficients cj cn agrees with the prediction made by Catani [15] 
using an infrared factorization formula. We use this formula to organize the 
finite part into terms arising from the expansion of the pole coefficients and 
a finite remainder, 

c (2),fi„ = c (l),ren _ ^ ^ren _ j (2 ) ^0) > (g) 

for i = {2,4,6,12}, and with the one- and two-loop insertion operators 
lW(e) and l( 2 )(e) given in [15]. 

As an example, we present our result for n/ A-contribution to the finite 

part c^' fin at two loops, 

«r*<*..*») = -M'i^t^ + 1 111(12)2 n 21 " 2 ?' 12 ' w 

V {X1+X2Y 4 Xi(l-X2) 

1 C(2) 1 13xi 2 +36xi-10xix 2 -18x 2 +31x 2 2 

+ 12(l-x 2 )xi ~ 18 (x 1 +x 2 ) 2 x 1 {l-x 2 ) n[ - X2 > 



xi-x 2 -2x 1 +4x 2r , , , 1 R(xi,x 2 ) 
H -, ■ T7 Ri(xi,x 2 N 



(xi+x 2 ) 4 ' 12xi(xi+x 2 ) 2 



5x 2 + 42xi + 5 



4 



(1+xO 2 l-3xi+3x! 2 „ xi 2 



4 — — - 72- 



1 — X 2 I — X1—X2 X1+X2 

l + 2xi 



+ 



1 1 6 



xi(xi+x 2 ) 2 



12xi(l-x 2 ) (xi+x 2 ) 2 

(Li.a-^-Li.a-x!))-^^-) 



-ll*n,N- H ^ 



2 1 xi(l-x 2 )' 
We have introduced the function R(xi, x 2 ), which is well known from [16], 
R(x!,x 2 ) = (10) 
Q ln(xi) ln(x 2 ) - ln(xi) ln(l-xi) + ^((2) - Li 2 (xi)) + (xi x 2 ) . 

In addition, it is convenient, to define the symmetric function Ri(xi,x 2 ), 
which contains a particular combination of multiple polylogarithms [11], 

Ri(xi,x 2 ) = (ln(xi)Lii 1 (— ^ — ,xi+x 2 ) - -£(2) ln(l-xi-x 2 ) (11) 
V \xi+x 2 J I 

+Li 3 (xi+x 2 ) - ln(xi)Li 2 (xi+x 2 ) - - ln(xi) ln(x 2 ) ln(l — x\ — x 2 ) 
-Lii 2 ( — v — ,x x +X2 ) - Li 2j i ( — ^ — ,xi+x 2 ) ) + (xi <-> x 2 ) . 



xi+x 2 / 

We have made the following checks on our result. As remarked, the in- 
frared poles agree with the structure predicted by Catani [15]. This provides 
a strong check of the complete pole structure of our result. In addition, we 
have tested various relations between the Cj. For instance, the combination 
x\c% is symmetric under exchange of x\ with x 2 . Finally, we could com- 
pare with the result for the squared matrix elements, i.e. the interference 
of the two-loop amplitude with the Born amplitude, and the interference of 
the one-loop amplitude with itself. The results of [5] are given in terms of 
one- and two-dimensional harmonic polylogarithms, which form a subset of 
the multiple polylogarithms [11]. Thus, we have performed the comparision 
analytically and we agree with the results of [5] . 



3. Conclusions 

Our result represents one contribution to the full next-to-next-to-leading 
order calculation of e + e~ — > 3-jets. It has been obtained by means of an 
efficient method based on nested sums and is expressed in terms of mul- 
tiple polylogarithms with simple arguments. As a consequence, our result 
can be continued analytically and applies also to (2 + l)-jet production in 
deep-inelastic scattering or to the production of a massive vector boson in 
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hadron-hadron collisions. At the same time, it provides an important cross 
check on the results for the squared matrix elements [5] with a completely 
independent method. 

After the results of section 2 had been presented at this conference, 
Garland et al. published results for the complete two-loop amplitude for 
e + e~ — > qqg. Our results are in agreement with ref. [17]. 
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